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Abstract 

The Kondo-Spin Glass competition is studied in a theoretical model 
of a Kondo lattice with an intra-site Kondo-type exchange interaction 
treated within the mean-field approximation, an inter-site quantum 
Ising exchange interaction with random couplings among localized 
spins and an additional transverse field T in the x-direction, which 
represents a simple quantum mechanism of spin flipping, in order 
to have a better description of the spin-glass state and in particular 
of the Quantum Critical Point(QCP). Taking here a parametrization 
T = aJ\ (where Jk is the antiferromagnetic Kondo coupling), we ob- 
tain two second order transition lines from the spin-glass state to the 
paramagnetic one and then to the Kondo-state. For a reasonable set 
of the different parameters, the two second order transition lines do 
not intersect and end in two distinct QCP. The existence of QCP in 
the Spin Glass-Kondo competition allows to give a better description 
of the phase diagrams of some Cerium and Uranium disordered alloys. 
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1. Introduction 



It is well known that there exists a strong competition between the Kondo 
effect on each site of a Kondo lattice and the magnetic ordering arising from 
the RKKY interaction between magnetic atoms in heavy fermion systems. 
The Doniach diagram[l] gives a good description of this competition: the 
Neel temperature T N is firstly increasing with an increasing of the antifer- 
romagnetic s-f exchange interaction constant Jr{> 0), then it is passing 
through a maximum and finally it tends to zero at the "quantum critical 
point" (QCP), with a second order transition at zero temperature. Such a 
decrease of T/v down to the QCP has been observed in many Cerium com- 
pounds, such &sCeAl 2 [2] or CeRh 2 Si 2 [3], under pressure. We also know that 
the Neel temperature starts from zero at a given pressure and incrases rapidly 
with pressure in Ytterbium compounds, such as YbCu 2 Si 2 [4\ or YbNi 2 Ge 2 [5], 
in good agreement with the Doniach diagram. Above the QCP, there exists a 
very strong heavy fermion character, but several possible behaviours, i.e. the 
classical Fermi liquid one with eventually a reduced Kondo temperature [6, 7] 
or different Non-Fermi-Liquid (NFL) ones, have been observed in Cerium or 
Ytterbium compounds [8, 9]. 

On the other hand, the disorder can yield a Spin Glass (SG) phase in ad- 
dition to the Kondo (mainly NFL) behaviour at low temperatures around 
the QCP in disordered Cerium or Uranium alloys. This is the case of the 
magnetic phase diagram of CeNii_ x Cu x alloys that has been extensively 
studied[10, 11], while the phase diagram of Ce 2 Aui^ x Co x Si 3 alloys presents 
the sequence of SG-AF-Kondo phases at low temperatures with increasing 
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Cobalt concentration[12]. The three phases, AF,SG, and NFL have been 
obtained at low temperatures for different concentrations in UCu 5 ^ x Pd x [13] 
or in U 1 - x La x Pd 2 Al 3 [14] disordered alloys. Thus, a SG-Kondo transition 
has been observed experimentally with increasing concentration in disordered 
alloys around the QCP. 

The purpose of our work here is to present a theoretical model that de- 
scribes the QCP for the spin glass-Kondo transition in Kondo lattices. The 
SG-Kondo transition was theoretically studied in a previous publication[15] 
and also in the presence of ferromagnetic ordering[16] or antiferromagnetic 
ordering[17], but the QCP was not described because we lacked a quantum 
mechanism that favored spin flipping. We present here an improvement of 
the previous model in order to obtain a good description of the QCP. 

More precisely, in previous publications[15, 16, 17] the resultant RKKY 
interaction was introduced by means of random, infinite range couplings 
among the S z components of the localized spins as in the Sherrington- 
Kirkpatrick (SK) model[18] and by neglecting the spin flip coupling of the 
transverse components. By using functional integral techniques with a static 
Ansatz in a replica symmetric theory, we obtained a magnetic phase dia- 
gram in the JkvsT plane, that showed the three different phases: paramag- 
netic, spin glass and Kondo. In spite of its complication, the model failed to 
describe a second order QCP at T = 0, because, by disregarding the spin 
flipping part of the Heisenberg Hamiltonian, we suppressed the tunneling 
mechanism, and magnetic ordering occurs only along the z-axis. In order 
to introduce a spin flipping mechanism and to avoid the intricacies of the 
random Heisenberg model[19], in the present paper the Heisenberg-like cou- 
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pling among the three spin components induced by the RKKY interaction is 
mimicked by a quantum Ising spin glass in a transverse field. It consists in 
an effective random interaction among the z-components, as we considered 
in[15], while the spin flipping part is simulated by a uniform transverse field 
in the x-direction. 

The infinite range quantum Ising spin glass in a transverse field T is one of 
the simplest models that presents a quantum critical point and it is equivalent 
to the model for a spin glass of quantum rotors[20, 21]. The dynamical 
properties at zero temperature are known and the existence of the QCP is well 
established [22], while the phase diagram in the TvsT plane has been studied 
by using the Trotter-Suzuki technique [23], and more recently [24] by the use 
of two fermionic representations of the spin operators in terms of Grassmann 
fields[25] that are more suitable to our purposes. As we did previously in 
reference [15], the Kondo effect will be studied in an approximation that is 
basically equivalent to the mean field decoupling scheme [6, 26]. The static 
Ansatz in the study of the spin glass transition is an approximation, similar 
to mean field theory, that is appropriate to describe the phase diagram. This 
is justified because in the model of M-components quantum rotors, that is 
exactly soluble for infinite M, the critical line is given by the singularity of the 
zero-frequency mode[20].The same occurs with the quantum Ising spin glass 
in a transverse field[22], where the singularity of the zero frequency mode 
determines the critical point. That is the reason why the static Ansatz, 
that describes only the zero frequency mode, describes the critical line in the 
phase diagram, although it would not give correctly the time dependence of 
the order parameter. 
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A related Hamiltonian has been considered in ref.[27] to describe NFL 
behaviour and a QCP in some heavy fermion compounds, although there are 
essential differences between this work and ours in the present paper and 
in ref.[15]. In ref.[27] it is considered a system of isolated Kondo impuri- 
ties, each one with a separate electron reservoir, and represented by the spin 
glass model of M-components quantum rotors in the limit of large M, when 
the problem is exactly soluble. The Kondo coupling provides the quantum 
mechanism and the transition line in the phase diagram is determined by 
the singularity of the zero frequency mode, displaying a QCP at zero tem- 
perature. The Kondo effect is described there [27] by isolated impurities and 
displays a continuous transition among different scaling regimes, in place of 
the sharp second order transition in Kondo latticesfl, 6, 26] 

This paper is organized as follows: in Sect. 2 we introduce the model, in 
Sect. 3 we discuss relevant results and we reserve Sect. 4 for discussions and 
conclusion. We refer the reader to ref. [15] for details in the mathematical 
calculations. 

2. The model 

We consider a Kondo lattice system with localized spins Si at sites % — 1 • • • JV, 
coupled to the electrons of the conduction band via a s-f exchange interac- 
tion. It is necessary to introduce explicitly the resultant RKKY interaction 
by means of a random, infinite range coupling among localized spins like in 
the Sherrington Kirkpatrick (SK) model for a spin glass [18]. To describe 
the Kondo effect in a mean-field-like theory it is sufficient to keep only the 
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spin-flip terms [6, 15, 26] in the exchange Hamiltonian, while the spin glass 
interaction is represented by the quantum Ising Hamiltonian and the trans- 
verse field T in the x-direction [24]. The transverse field T represents a simple 
quantum mechanism of spin flipping and mimics the more complicated trans- 
verse part of the Heisenberg Hamiltonian [19]. 
The Hamiltonian of the model is: 

H = H K + H SG (1) 

n K = J2 ^n ks + e J2 n{ + J K Y\ S Pli + S Ji4i\ (2) 

7-t-SG = - J K S fi S h ~ 2F 51 S fi ( 3 ) 
where J K > 0, S ± = S x ± iS y and 

Sfi = ^2 fis&s,s>fis> 
s,s' 

Sdi = d ls^s,s'dis' (4) 
s,s' 

f} s , fis{di S ,di) are creation and destruction fermion operators for electrons 
with s =t or | in the localized (conduction) band. We indicate by g_ the 
Pauli matrices and we have n ka = d\. a dk a the conduction electrons occupation 
number. The energies eo(efc) are referred to the chemical potentials 
respectively. 
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The coupling in Eq. (3) is an independent random variable [18] with 
gaussian distribution of zero mean and variance < J?- >= 8J 2 /N . 
Functional integration techniques have proved to be a suitable approach to 
describe phase transitions in disordered quantum mechanical many-particle 
systems [25]. The static approximation within this formulation consists in 
neglecting time fluctuations of the order parameter, and when it is combined 
with the neglect of space fluctuations it leads to the usual Hartree-Fock, mean 
field like approximation. When dealing with the Hamiltonian in Eq. (1)- 
Eq.(3), we notice that in the limiting case Jk = we obtain a quantum Ising 
spin glass in a transverse field that has been studied with the static approxi- 
mation in [24], while for J = we recover the mean field approximation that 
has been used successfully to describe the Kondo lattice [6, 8, 26]. We follow 
closely the formalism of Ref. ([15]) to write a Lagrangian formulation [25] in 
terms of anticommuting, complex Grassmann variables ipi S (cu) , ipi S (cu) asso- 
ciated to the conduction and localized electron fields, respectively, together 
with the spinors 

v <fil(u) ) 

where u = (2n + 1)tt are the Matsubara frequencies. The partition function 
is now expressed 

Z = J Z%V)^V)e A (6) 

We now follow the same steps as in Ref. ([15, 16, 17]) and in the static, mean 
field like approximation the action A may be written 

A = A + A K + A SG (7) 



(5) 



with Aq being the action for non-interacting electrons in a transverse mag- 
netic field: 

A o = J2J2 ^[M [(^ - fco + + tplfaSij - Pti^tpAu)] (8) 

while the Kondo part of the action is decoupled in the mean field approxi- 
mation as in Ref. ([15]) 

a k = ^EE4(^)^H][E^U-H] (9) 

Here (3 is the inverse absolute temperature and we also have 

A 8a = T l J » s fi s h ( 10 ) 

where in the static approximation [15, 16, 17, 24] 

The Kondo order is described by the complex order parameter 

A. = -^E<rf>M.M> (i 2 ) 

i,u> 

that in a mean field theory [6, 15, 26] describes the correlations X\ = (fj s di S ) 
and A s = (d\ s f is ). The approximation used in Eq.(12), which is equivalent 
to the "slave boson" method, is certainly one of the best practical methods 
used for the Kondo lattice problem. The Kondo temperature for the lattice is 
determined here by the temperature at which A becomes equal to zero and 
this approximation gives a fairly reasonable description of the Kondo phase. 
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However, it is unable to give a "mixed" SG- Kondo phase, as it was previously 
shown in the case of the antiferromagnetic-Kondo transition [28]. Except 
for the presence of the transverse field T in A in Eq. (8), the calculation 
follows the same steps as in Ref. ([15]) and we refer the reader to this paper 
for details. We show in Ref. ([15]) that standard manipulations give for the 
averaged free energy within a replica symmetric theory: 

(3F = 2(3 J K \ 2 + + 2qx) - Pn (13) 



where 



1 f+oo y n r+oo ^ 

Pn = H^ U Dz iU n^exp(^ln|G-.»|)-l} (14) 



and q,% and A must be taken at their saddle point value. Here a is the 
replica index and q is the SG order parameter[15, 18, 24] while x — fix is the 
static uniform spin susceptibility of the localized f-electrons. 
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9a« = ? + X = l™-E< 5 ^> > ( 16 ) 

n ^ a 

j 1 2 

We use the notation Di = -j§^e~^ x and we refer the reader to Ref. ([15]) 
for the details of the calculation. The matrix G ija (uj) in Eq. (14) is the time 
Fourier transform of the matrix Green's function GV,-(t) = i{Ttpi(T)ip^(0)) for 
the localized electrons in the presence of random fields Zj and £ Q j at every 
site, and satisfies the equation[15] 

G^Klu) = [iw - (3e - aji ja + (3To^}5tj - (5 2 JlAU^Hl (17) 



where 




while 7ij(a>) is the time Fourier transform of the conduction electron Green's 
function 7y(r) = i(Tip is (T)<fi S (0)) and is given by 

1^ = N - M&a - PUj (19) 

We obtained in Eq. (17) the Green's function for the f-electrons in a Kondo 
lattice[15], but now in the presence of a random field hj at every site that 
prevents us from diagonalizing by means of a Fourier transformation. In the 
pure SG limit Jk = 0, the Green's function in Eq. (17) is local and the 
integrals in Eq. (14) reduce to a one site problem, while in the Kondo limit 
J = the random fields vanish and the integrals separate in reciprocal space, 
giving the known results for a Kondo lattice. We adopt here a decoupling 
approximation introduced in Ref. ([15]) that reproduces correctly these two 
limits, and corresponds to consider independent Kondo lattices, in place of 
the independent Kondo impurities in Ref. ([27]). We replace the Green's 
function G_ i j a {uj ) hi a ...hj a ...h Na ) by the Green's functions G^ u (ui,hj a ) with, 
j = 1...JV, of iV independent Kondo lattices, each one with a "uniform" field 
hj a at every site /i, v by means of the approximation 

j 

where G F (w,/ij a ) is the f-electron Green's function for a fictitious Kondo 
lattice that has an uniform field hj a at every site /i, v and satisfies the equa- 
tion: 

G£(u, h ja ) = [{iu - Pe )l - a z h ja + pTeJS^ - f3 2 J|A Vl (21) 
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where from Eq. (19) 



Now Eq. (21) may be easily solved by a Fourier transformation with the 
result: 

In |GJ>, h ja ) | = £ln l<5>, M I (23) 

where 

G^(o;, fc ja ) = [iw - /3e ]l - ^/i jQ + f3Fa x - f3 2 J 2 K \ 2 . 1 - . (24) 
— — *o; - pe fc 

We may now introduce Eq. (24) and Eq. (20) in Eq. (14), the integrals over 

the fields separate and we obtain 

/+oo c+oo 1 

Dz\n{ i^exp(- (25) 
-oo J— oq iV _ _ 

fc * 

with 

£S.(fc,lf) = X;M<V 1 (",'0] (26) 
and /i is given in Eq. (18), with z and £ in place of Zj and £j Q , while 

H = yV + {pry (27) 

The sum over the fermion frequencies is performed in the standard way by 
integrating in the complex plane [25], with the result: 

S,(k, H) = In [(1 + e~"s+ ) (1 + e~ u >- )] (28) 

where 

uj s ± = l -[f3e k + sH] ± {- A {l3e k - sH) 2 + (f3J K \) 2 } 1 *. (29) 
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We consider eo = that corresponds to an average occupation (rif) = 1, per 
site. Replacing sums by integrals , in the approximation of a constant density 
of states for the conduction band electrons, p(e) = p = for — D < e < D, 
we obtain from Eq. (26) to Eq. (29) the final expression for the free energy 
in Eq. (13): 

1 r+oo r+oo 

f3F = 2(3 J K \ 2 + -(3 2 J 2 (x 2 + 2qx) ~ Dz\n{ D^} (30) 

Z J -co J— oo 

with 

1 r+/3D ( r _L If) 

E(H) = 4t: dx In {cosh +CO sh(v / A)}, (31) 

(3D J-/3D 2 

^ = \{x-H) 2 + {(3J K \f (32) 
and from Eq. (27) we have 

H = ^((3ry + ((3jy(^2qz + ^) 2 (33) 

3. Results 

The SG-Kondo transition is described in the present model by the three 
parameters Jk, J and T, which cannot be considered here as completely 
independent from each other. In fact, the Kondo effect and the RKKY in- 
teraction originate from the same intra-site exchange interaction, but the 
necessity of considering an additional inter-site exchange term has been al- 
ready recognized [6, 28] and a relationship giving the intersite exchange pa- 
rameter varying as (Jr) 2 has been introduced in ref.[6]. Since the transverse 
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field is introduced here to mimic the spin flipping part of the Heisenberg 
Hamiltonian that originate in the RKKY interaction, we assume a similar 
relationship T = a(J K ) 2 , a < 1 to have a better description of the SG-Kondo 
transition. In particular, the consideration of the previous relationship for 
T will be important for the existence of the QCP and the comparison with 
experiment. The saddle point equations for x, q and A are obtained by ex- 
tremizing f3F with respect to these order parameters. By making q = A = 
in the saddle point equations we obtain the equations for the second order 
critical lines Ti(J K ) and T 2 (Jk) that separate the paramagnetic phase from 
the spin-glass and Kondo phases, respectively. They are given by: 

I r+oo 

*(® = ?— T / D ^ ^HH )/H (34) 

1 + Jn J— 00 



where 



/+00 
D£ cosh(# ) (35) 
-00 

and H = H(q = 0). 

Thus we obtain for the spin glass transition temperature Ti(J K )[2A\. 

^ - (36) 

that together with Eq. (34) gives Ti and x(/3i), while the Kondo transition 
temperature T 2 (Jk) is given by 



k 



+00 



4(l + J )^^ COSh W 2) 

xJLr*°«fa ±_^Sll = (37) 
(3 2 DJ-p 2 D cosh(x/2) ^£ v > 
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that must be solved together with Eq. (34). 

The numerical solution of the saddle point equations as a function of j and 
^f- yields the phase diagram in Fig.l, where we considered D/J = 12 and 
two values of aJ. The two cases corresponding to these two values of aJ 
are given in Fig.l . The first case (represented by full lines) corresponds to 
the really interesting situation of the phase diagram with two QCP, while the 
second case (represented by the dashed line T') is plotted here for comparison 
to show the case without QCP. 

In the first case, corresponding to aJ = 0.01348 (solid line) the second order 
critical line Ti(J^) that separates the paramagnetic from the spin glass phase 
ends at a QCP corresponding to a Jk value called J^x here, while the second 
order critical line T 2 (Jk) that separates the paramagnetic from the Kondo 
phase ends at a second QCP J C K2 an d we have here J^ 2 > Jfci- 
An asymptotic calculation at low temperatures gives: 

k B T x = ^-[2y/2J - aJ 2 K ] (38) 

and the corresponding value for J C K1 at the QCP is given by a( J^i) 2 — 2\[23 . 
On the other hand, close to the second QCP an asymptotic calculation gives: 

^^^' hS/j ' (39) 

where T K = D exp (— 2D/ J K ). Thus the value of J£- 2 is the solution of 
T k (Jk 2 ) = ®(Jk2) 2 i n the asymptotic limit of very large D values. In order 
to obtain a solution J^ 2 that makes T 2 = in Eq. (39) it is necessary that 
AaD = 0.6473, what gives aJ = 0.013485 for D = 12 J, and for these values of 
the parameters in Fig.l, we have obtained two QCP with J^ 2 > J^ 1 . In the 
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second case corresponding to a smaller value of a, a J = 0.1344, the two QCP 
disappear and we obtain a single transition line ( dashed line T') separating 
the paramagnetic from the ordered ( spin-glass and Kondo) phases. The 
physical case shown in Fig.l, where the two QCP are very close to each other 
and almost equal, corresponds well to experimental phase diagrams obtained 
in some Cerium or Uranium disordered alloys. We obtain from Eq. (34) in 
the disordered region that x — > 0.5 when j3 — * 0, thus giving the expected 
Curie- Weiss behaviour for the static susceptibility x — Px, while for (3 — > oo 
we get x —> r + Vr 2 - 8 J 2 when r = aJ\ > 2y/2J . We recover here the 
square root singularity in the linear static susceptibility that is present in 
other quantum spin glass models[20, 21, 22]. It is interesting to remark that 
X = xT = S(S + 1) is the residual effective localized spin in Curie- Weiss 
theory. The value for x on the critical line T 2 (Jk) is presented in Fig.l, 
where we notice the steep drop to almost zero close to the QCP J C K2 . 

4. Conclusions 

We study in this paper the phase transitions in a heavy fermion system rep- 
resented by a Hamiltonian that couples the S z components of the localized 
spins of a Kondo lattice [7] with random, long range interactions, like in the SK 
model for a spin glass[18], while the transverse components are acted upon 
by a field in the x-direction[24]. As the transverse field T mimics the spin flip 
part of the Heisenberg coupling among localized spins, that originates in the 
RKKY interaction, we assume that r rs olJ\ where Jk is the antiferro mag- 
netic Kondo coupling[6]. Using functional integrals techniques and a static, 
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replica symmetric Ansatz for the Kondo and spin glass order parameters, 
we derive a mean field expression for the free energy and the saddle point 
equations for the order parameters. The use of the static Ansatz in the case 
of the transverse spin glass is justified, because it is the singularity of the 
zero frequency mode that determines the critical line[20, 21, 22]. Numerical 
solution of the saddle-point equations allows us to draw the magnetic phase 
diagram in the JkvsT plane ( presented in Fig.l) for fixed values of J and D 
with D/J =12 and for two values of aJ, aJ = 0.01348 or aJ = 0.01344. 
Figure 1 shows three distinct phases. At high temperatures, the "normal" 
phase is paramagnetic with vanishing Kondo and spin glass order parame- 
ters, i.e. A = q = 0. When temperature is lowered, for not too large values 
of the ratio Jk/ J, a second-order transition line is found at T = T\ to a spin 
glass phase with q > and A = O.The critical line Ti(J K ) ends at a quantum 
critical point J K1 . Finally, for large values of the ratio Jr/ J > JkiI J ■> where 
Jr2 — Jrii there is a second order transition line to a Kondo state with 
A > 0. These results are very sensitive to the value of D/J and a J. For 
D/J = 12 and a J < 0.01344 the QCP disappear while they are favored for 
larger values of a J. We can also remark that we get here only "pure" Kondo 
or SG phases and never a mixed SG-Kondo phase with the two order pa- 
rameters different from zero as already noticed for the competition between 
Kondo and antiferromagnetic phases [28]; this result is probably connected 
to the approximations used here to treat the starting Hamiltonian. A QCP 
has been observed experimentally in several Cerium and Uranium disordered 
alloys[12, 13, 14] and the phase diagram shown in Fig.l which yields two QCP 
improves considerably the description of the spin glass-Kondo transition with 
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respect to previous publications[15, 16, 17]. However, the experimental sit- 
uation is generally not really clear; for example, there is no experimental in- 
formation on the precise nature of the SG-Kondo transition in CeNii_ x Cu x 
alloys [10]. Moreover the phase diagrams of several systems involve also an 
antiferromagnetic phase and this case is theoretically studied elsewhere[17]. 
Our theoretical results describe also the spin glass and Kondo phases in Ura- 
nium alloys. An unsolved basic question concerns also the existence or not 
of a "mixed" SG-Kondo phase in Cerium and Uranium disordered alloys and 
this problem is worth of being studied experimentally in more detail. Thus, 
further experimental work is necessary, but our model yields a new striking 
point in the behaviour of heavy fermion disordered alloys in the vicinity of 
the quantum critical point. 
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Figure captions 



Figure l:Phase diagram in the T — Jk plane as a function of T / J and Jk/ J 
for the transverse field T = aJ K , a J = 0.01348 ( solid line) and a J = 0.01344 
(dashed line), for the conduction bandwidth D/ J = 12. The critical (solid) 
second order line Ti for low values of Jk separates the paramagnetic phase 
(q — X — 0) for high temperatures from the spin-glass phase (q > 0, A = 0) 
at low temperatures and ends at a QCP. The critical second order (solid) line 
T 2 for large values of Jk separates the paramagnetic phase from the Kondo 
phase (A > 0, q = 0) and ends at a second QCP. The critical ( dashed) line 
T'does not have a QCP. The dash-dotted line represents the "pure" Kondo 
temperature Tk and the dotted line represents the residual magnetic moment 
X c {Jk) on the line T = T 2 . 
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